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Abstract. Given a strictly hyperbolic, genuinely nonlinear system of conservation laws, we prove 
the a priori bound \\u(t, •) — u E (t, -)|| L i = + t) ■ y/e\ lne| on the distance between an exact 

BV solution u and a viscous approximation u £ , letting the viscosity coefficient e — > 0. In the proof, 
starting from u we construct an approximation of the viscous solution u e by taking a mollification 
u * (p and inserting viscous shock profiles at the locations of finitely many large shocks, for 
each fixed e. Error estimates are then obtained by introducing new Lyapunov functionals which 
control shock interactions, interactions between waves of different families and by using sharp decay 
estimates for positive nonlinear waves. 



1 - Introduction 

Consider a strictly hyperbolic system of conservation laws 



Here A(u) = Df(u) is the Jacobian matrix of /. Given an initial data u(0,x) = u(x) having small 
total variation, the recent analysis in [BiB] has shown that the corresponding solutions u £ of (1.2) 
exist for all t > 0, have uniformly small total variation and converge to a unique solution of (1.1) 
as e — > 0. The aim of the present paper is to estimate the distance ||u £ (i) — u(t) || L1 , thus providing 
a convergence rate for these vanishing viscosity approximations. 





together with the viscous approximations 



u\ + A(u £ )u% = eu% 



(1.2) 



We use the Landau notation 0(1) to denote a quantity whose absolute value remains uniformly 



bounded, while o(l) indicates a quantity that approaches zero as e — > 0. Our main result is the 
following. 

Theorem 1. Let the system (1.1) be strictly hyperbolic and assume that each characteristic field is 
genuinely nonlinear. Then, given any initial data u(0, •) = u with small total variation, for every 
t > the corresponding solutions u,u £ of (1.1) and (1.2) satisfy the estimate 

\\u £ {t, •) - u(t, •) || l1 = 0(1) • (1 + r)Ve| lne| Tot.Var.{u} . (1.3) 



Remark 1. For a fixed time r > 0, a similar convergence rate was proved in [BM] for approximate 
solutions generated by the Glimm scheme, namely 

||« GHmm (r,-)-n(r,-)|| Ll =o(l)- V£|lne|. 
Here e « Ax rj At measures the mesh of the grid. 

Remark 2. For a scalar conservation law, the method of Kuznetsov [K] shows that the convergence 
rate in (1.3) is 0(1) • e 1 / 2 . As shown in [TT], this rate is sharp in the general case. 

In the case of hyperbolic systems, in [GX] Goodman and Xin have studied the viscous approxi- 
mation of piecewise smooth solutions having a finite number of non-interacting shocks. With these 
regularity assumptions, they obtain the convergence rate 0(1) • e 7 for any 7 < 1. On the other 
hand, the estimate (1.3) applies to a general BV solution, possibly with a countable everywhere 
dense set of shocks. 



To appreciate the estimate in (1.3), call St and S £ the semigroups generated by the systems 
(1.1) and (1.2) respectively. As proved in [BCP], [BLY] and [BB], they are Lipschitz continuous 
w.r.t. the initial data, namely 

\\S t u - Stv^ < L \\u - u|| Ll , (1.4) 

||^-^|| L1 <L||«-i;|| Ll . (1.5) 
The Lipschitz constant L here does not depend on t,e. By (1.4), a trivial error estimate is 

\u £ (t + h)- S h u £ <+^ 1 



||« e (r)-«(r)|| Ll = L-J Q | fe hm + l|M + - ^1^1 j dt = L . jf \\e U £ xx (t)\\^ dt . 

However, ll^^Ci) || L i grows like e _1 , hence the right hand side in the above estimate does not 
converge to zero as e — > 0. 



2 



We thus need to take a different approach, relying on (1.5). Let e > be given. It is well 
known (see [B2]) that one can construct an e'-approximate front tracking solution u of (1.1), with 

\\u(0) — u|| Ll < e', \\u(t) — u(r)|| L1 < e', 

and such that the total strength of all non-physical fronts is < e'. Here we can take for example 
e' = e -1 / £ . Since the errors due to the front tracking approximation are of order e' « e, in the 
following computations we shall neglect terms of order 0(1) • e 1 as they can be made arbitrarily 
small by a suitable choice of e' . For sake of definitiness, we shall always work with the right- 
continuous version of a BV function. Since all characteristic fields are genuinely nonlinear, it is 
convenient to measure the (signed) strength of an i-rarefaction or of an i-shock front connecting 
the states u~ , u + as 

a = Xi(u + ) - Xi(u~) , 

where Aj denotes the i-th eigenvalue of the matrix A(u). We follow here the notations in [B2], and 
call 

v ( u ) = ^2Wa\, Q(u)= (1-6) 

respectively the total strength of waves and the interaction potential in a front tracking solution u. 
The second summation here ranges over the set A of all couples of approaching wave fronts. 

For notational convenience, we shall simply call u the e'-approximate front tracking approx- 
imation, also assume that u = u(0) is piecewise constant. Since e' « e, this will not have any 
consequence for our estimates. In the sequel, we shall construct a further approximation v = v(t, x) 
having the following properties. 

Let = to < ^i < • • • < In = t be the interaction times in the front tracking solution u. Then 
v is smooth on each strip [U-i , U[ xM. Moreover, calling Sq = Tot.Var.ju}, one has 

\\v(0) - u\\^ = 0(1) ■ SoVi , ||i;(t)-u(t)|| l1 =0(1)-W£, (1-7) 

j \v t + A(v)v x -ev xx \dxdt = O(l)-8 (l + T)^\\ns\, (1.8) 

I \v{t i ,x)-v{t i -,x)\dx = O{l)-5 ^\\ne\. (1.9) 

l<i<N J 

Having achieved this step, by the Lipschitz continuity of the semigroup S\ in (1.5) we can then 
conclude 

IK(r) - «(r)|| Ll < ||^« - «(r)|| Ll + \\v(r) - u(r)\\^ 

< L \\u — v (0)|| Ll + L / / \vt + A(v)v x — ev xx | dxdt 

(1-10) 

+ L 2J ll W (^' X ) ~ V ( t i-^ X )\\l J i + |K T ) ~ n ( T )|| L 1 

l<i<N 

= C(l)-5o(l + r)v / i| hie|. 



To construct the approximate solution v, we first consider a mollification of u w.r.t. the space 
variable x. Let ip : M <— > [0, 1] be a smooth function such that 



<p{s) = 



if \s\ > 3 



s<p'(s)<0, ip(s) = ip(-s), 



J (p(s) ds = 1 



For 5 > small, define the rescalings <ps(s) = 5 1 (p(x/5) and the mollified solutions v 5 (t) 
u(t) * tps, so that 

v s (t, x) = j u(t, y)(fs(x -y)dy. 
Recalling that 5 = Tot.Var.j'u}, one has 



for all t>0. 



(1.11) 



dx 



(1.12) 



(1.13) 



Tot.Var.{u(i)}, ||<(i)|| Ll = 0(1) ■ 8 , 
We now observe that 

IK ~ = J J { u ( x )- u (y)) L Ps(x-y)dy 

< J Tot.Var.ju; [x - 5, x + 5]} dx = O(l)-S S. 
To estimate the distance between v 5 and u £ , we first compute 

y | £U L( 2; )| dx = e j \( u x * (P5,x)(x)\dx < £||« x || L i • ||<A5,z|| L i = O(l) • #0 | 

J \v s t +A(v 5 )v s x \dx = J J (A(v s (x))u x (y)~A(u(y))u x ( y )^(x-y)dy 
< J ^ \a(v s (x)) - A(u(y))\<p(x - y) dx^j \u x {y)\dy 

= 0(1) ■ \\DA\\ C o j Osc.ju; [y - 6, y + 5}} \u x (y)\ dy . 

For simplicity, the formulas (1.13)-(1.14) are here written in the case where the function u is 
absolutely continuous. In the general case, the same estimates hold, by replacing \u x \dx with the 
measure \D x u\ of total variation of u G BV . 

If u is a Lipschitz continuous solution of (1.1), the oscillation of u on any interval of length 25 
is 0(1) • S. Hence, performing the above mollifications, we would obtain 



dx 



(1.14) 



/ 



(1.15) 



\v* + A(v 5 )v 5 x \dx = O(l)-8 6 . 

Choosing 5 = y/e, by (1.12)— (1.15) we thus conclude 
IK(r) - u(t)|| l1 < \\S*u - t/(r)L + ||«*(r) - u(t)|| l1 

<L\\u-v 5 (0)\\^+L J \ v f + A(v 5 )v 5 x -ev 5 xx \dxdt+\\v 5 (T)-u(r)\\ hl 
= O(1)-5 (1 + t)VI. 



(1.16) 



In general, however, the solution u is not Lipschitz continuous. The best one can say is that u is 
a function with bounded variation, possibly with countably many shocks. Hence the easy estimate 
(1.16) does not hold. For genuinely nonlinear systems, the additional error terms due to centered 
rarefaction waves can be controlled by carefully estimating the decay rate of these waves. Error 
terms due to small shocks will be estimated by suitable Lyapunov functionals. However, there is 
one type of wave-fronts which is responsible for large errors in (1.14), namely the large shocks of 
strength » \fe. In a neighborhood of each one of these shocks, a more careful approximation is 
needed. Instead of a mollification, we shall insert an approximate viscous shock profile. 

Our construction goes as follows. By the same argument as in [BC1] (see Proposition 2 on 
p. 17), given p > one can select a finitely many shock fronts 

t^x a (t) t € T a =]t~, t+[, a = l,...,u, 

with v = 0(1) • So/p, having the following properties. 



+ 



• For every t £T a (apart from finitely many interaction points) the left and right states u~ , u 
are connected by a shock, say of the family k a , with strength |cr a (i)| > p/2, while |u Q (f*)| > p 
for some t* G T a . Moreover, every shock in the front tracking solution u with strength > p is 
included in one of the above fronts. 

For each a and each t 6 T a (apart from finitely many interaction points), let co a be the viscous 
shock profile connecting the states u~,u+. Calling X a the shock speed, we thus have 

u'a = (A(u a ) - \ a )uj' a , lim u a (s) = . 

We choose the parameter s so that the value s = corresponds roughly to the center of the 
travelling profile. This can be achieved by requiring 

/0 poo 
\uo a (s) - u~\ ds = / \u a (s) - ds . (1-17) 
-oo J 

For the system (1.2) with e-viscosity, the corresponding rescaled shock profile is s i— > 0J £ a (s) = 
oj a (s/e). On the open interval 

Jocit) = ]x a (t)-S , x a (t) + S[ 
we now replace the mollified solution by a shock profile. Define the functions g a , £o a , by setting 

Vs(y)dy + u- / <ps(y)dy, (1.18) 

-oo Jf 



f c^(<XO) if 5[, 

<_ if £>5, (1-19) 



where 



re 



4(^-0 



if iei<^, 
if ~^<e<-^ 



(1.20) 



4(v^+5) 

Notice that w Q is essentially an e-viscous shock profile, up to a C 1 tranformation that squeezes the 
whole real line onto the interval J a (t). Moreover, g a is the mollification of the piecewise constant 
function taking values u~, u+ with a single jump at x a . The above definitions imply that Cj a = g a 
outside the interval J a {t). Finally, for every i>0we define 



V = U * ip S + ^ foot - Qa), 



X21) 



a€BS 



where the summation ranges over all big shock fronts. In the remainder of the paper we will show 
that, by choosing 

S = y/e, p = 4y/e\lne\, (1.22) 



all the estimates in (1.7)-(1.9) hold. By (1.10), this will achieve a proof of Theorem 1. 



2 - Estimates on rarefaction waves 

Throughout the following we denote by Ai(«) < • • • < \ n (u) the eigenvalues of the A(u) = 
Df(u). Moreover, we shall use bases of left and right eigenvectors k(u), r^(u) normalized so that 

VA^-nM^l, ii(«)-r i («) = {j !| ly.] (2.1) 

According to (1.14), outside the large shocks we have to estimate the quantity 

E{t) = £ J Osc.{u; [y-S, y + 5}} \u x (y)\ dydt . (2.2) 

Centered rarefaction waves can have large gradients, and hence give a large contribution to the 
above integral. However, for genuinely nonlinear families, the density of these waves decays rapidly, 
as t -1 . We now give an example where the integral (2.2) can be easily estimated. 

Example 1. Assume that the solution u consists of a single centered rarefaction wave of the i-th 
family (fig. 1), connecting the states u~ ,u + . Call s u(s) the parametrized i-rarefaction curve, 
so that 

d) = ri(u>), oj(0) = u~ , w(c) = u + 
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for some wave strength a > 0. We then have 



u 



u~ if x/t < Xi(u~) , 

(t,x) = ^ u(s) if x/t = Aj(w(s)) sG[0, cr], 

■u+ if x/t > Xi{u + ). 



If X is an upper bound for the length of all eigenvectors ^(u), we have 

Osc. {u(t); [y-S, y + 5}} < K-mm{a, 25/t] , j \u x (x)\dx < Ka . 

Hence the quantity in (2.2) satisfies 

E(t) < / K 2 a 2 dt + / K 2 a—dt = 2K 2 5a ( 1 + In ) • (2-3) 
The choice (5 = y/i would thus give the correct order of magnitude 0(1) ■ \fe \ \ne\ ■ Tot.Var.{«}. 





figure 1 



figure 2 



Of course, a general BV solution of the system of conservation laws (1.1) is far more complex 
than a single rarefaction. It can contain several centered rarefactions originating at t = and also 
at later times, as a result of shock interactions (fig. 2). Moreover, the crossing of wave fronts of 
other families may slow down the decay of positive waves. Nevertheless, the forthcoming analysis 
will show that, in some sense, Example 1 represents the worst possible case. Using the sharp decay 
estimate for positive waves in [BY] and a comparison argument, we shall prove that the total error 
due to steep rarefaction waves for an arbitrary weak solution is no greater than the error computed 
at (2.3) for a solution containing only one centered rarefaction. In the present section, all the 
analysis refers to an exact solution. A similar result can then be easily derived for a sufficiently 
accurate front tracking approximation. 

We begin by recalling the main results in [BY]. Given a function u : M n with small total 

variation, following [BC] and [B2], one can define the measures // of i- waves in u as follows. Since 



7 



u 6 BV, its distributional derivative D x u is a Radon measure. We define fi % as the measure such 
that 

^ = k(u) ■ D x u (2.4) 
restricted to the set where u is continuous, while, at each point x where u has a jump, we define 

^({x})=a l , (2.5) 

where <jj is the strength of the z-wave in the solution of the Riemann problem with data u~ = u(x—), 
u + = u(x+). In accordance with (2.1), if the solution of the Riemann problem contains the 
intermediate states u~ = u , lo\, . . . , uj n = u + , the strength of the i-w&ve is defined as 

<7j = Xi(cOi) - Aj(a>j_i). (2.6) 

Together with the measures \i l we also define the Glimm functionals 

v( U ) = j2\A(m, 

i 

Q( u ) = (1^1 l^Difav)' x <y} (^l^Dlix^); x^y}, 

i<j i 

measuring respectively the total strength of waves and the interaction potential. 
We call respectively the positive and negative parts of //, so that 

^ = / u l+ -//", |//| =//+ + //". (2.7) 

In [BY], the authors introduced a partial ordering within the family of positive Radon measures: 

Definition 1. Let be two positive Radon measures. We say that y < y! if and only if 

sup n(A) < sup y'{B) for every s > . (2.8) 

meas(A)<s meas(B)<s 



Here meas(j4) denotes the Lebesgue measure of a set A. In some sense, the above relation 
means that // is more singular than y. Namely, it has a greater total mass, concentrated on regions 
with higher density. Notice that the usual order relation 

y < // if and only if v(A) < n'{A) for every A C M 

is much stronger. Of course fi < // implies /j, ^ //, but the converse does not hold. 

Given a solution u of (1.1), we denote by y\ + the measure of positive i-waves in u(t, ■). In 
particular, ^ refers to the positive z-waves in u at the initial time t = 0. An accurate estimate 
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of these measures is obtained by a comparison with a solution of Burgers' equation with source 
terms. 



Proposition 1. For some constant k > and for every small BV solution u = u(t,x) of the 

system (1.1) the following holds. Let w = w(t,x) be the solution of the Cauchy problem for 
Burgers' equation with impulsive source term 

w t + (w 2 /2) x = -Ksgn(x) ■ j t Q(u(t)) , (2.9) 

w(0,x) =sgn(x) ■ sup . (2.10) 

meas(A)<2\x\ * 

Then, for every t > 0, 

ti + <D x w{t). (2-11) 



For a proof, see [BY]. 

The ordering relation (2.8) can be better appreciated in terms of rearrangements. More 
precisely, let /ibea positive Radon measure on M, so that \i = D x v is the distributional derivative 
of some bounded, non-decreasing function v : 1R i— > M. We can decompose 



as the sum of a singular and an absolutely continuous part, w.r.t. Lebesgue measure. The absolutely 
continuous part corresponds to the usual derivative z = v x , which is a non- negative L 1 function 
defined at a.e. point. We shall denote by z the symmetric rearrangement of z, i.e. the unique even 
function such that 

z(x) = z(—x) , z(x) > z(x') if < x < x 1 , 

meas^jx; z(x) > c}^ = meas^jx ; z(x)>c}^j, for every c > . 

Moreover, we define the odd rearrangement of v as the unique function v such that 

v(- x ) = -v( x ) , {}(o+) = ^ sing (iR) , 



r 

v(x) = v(0+) + / z{y) dy for x > . 
Jo 



'o 

By construction, the function v is convex for x < and concave for x > 0. We now have 

Proposition 2. Let = D x v and y! = D x v' be positive Radon measures. Call v,v' the odd 
rearrangements ofv,v', respectively. Then /i ^ D x v ^ y. Moreover 

v(x) = sgn(x) • sup . (2.12) 

meas(A)<2\x\ 2 
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Moreover, 



H ^ p' */ and only if v{x) < v'(x) for all x > . 



(2.13) 



The relevance of the above concepts toward an estimate of the quantity in (2.2) is due to the 
next three comparison lemmas. 

Lemma 1. Let u : M i— > M be a non- decreasing BV function and let u be its odd rearrangement. 
Then 

/CO n CO 

Tot.Var.ju; [x - p, x + p]} du(x) < 3 / [u(x + p) - u(x - p)] du(x) . (2.14) 
-OO J — CO 



Proof. We begin by defining a measurable map x i— > <p{x) from M onto 1R+ with the following 
properties. 

(i) ip(x) = for all points x in the support of singular part of the measure u x . 

(ii) u x (x) = ii x (ip(x)) for every x where u is different iable. 
(Hi) meas(<p _1 (^4)) = 2meas(^4) for every A C M + . 

We now have 

/CO 
Tot.Var.ju ; [x — p, x + p]} du(x) 
■ CO 

= I / + / + p) — u(x + p)] (fu(s 

\J<p{x)<p Jip(x)>pJ 

= h + h- 

We now estimate I\ and /2 separately as follows. 



[X 



h = [u(x + p) — ii(x — p)] d-u(x) 

J<p{x)<p 

< / 2u(p)du(x) 

J<p{x)<p 

<A{u{p)) 2 
f P 

< 2 / [u(x + p) — ti(x — p)] du{x). 

J-p 



(2.15) 
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I<2. < / / [u x (x) D x u(x + s)] dsdx 

■J <p(x)>p J —p 

/■OO 

< 4p / [u x (x) D x u(x - p)] dx 

J P 

POO 

= 4p/ ^(x + p)^(x) (2.16) 
./o 

< 2 / [u(x + p) — n(x — p)] dii(x) 

Jo 

/OO 
[u(x + p) — -u(a; — p)] du{x). 
-OO 

For x > p, we are here using the inequality 

2pu x (x) < u{x) — u(x — 2p). 

Moreover, calling /, g the non-increasing even rearrangements of two positive, integrable functions 
/, g, one always has 

/OO rOO 
f{x)g{x)dx< / f(x)g(x)dx. (2.17) 
-OO J — OO 

Together, (2.15) and (2.16) yield (2.14). □ 

Lemma 2. Let v,w be two non- decreasing BV functions. If D x v ^ D x w then the odd rearrange- 
ments v, w satisfy 

/OO 
[w(x + p) — w(x — p)] dw(x) . (2.18) 
-OO 

Proof. By an approximation argument, we can assume that v and w are smooth. Without loss 
of generality, we can assume 0(±oo) = i/)(±oo). By assumptions, v(x) < w(x) for all x > 0. We 
consider a parabolic equation with smooth coefficients 

z t = a(t,x)z xx , (2.19) 

with a(t, x) = a(t, —x) > 0, having a solution such that 

z(0,x) = w(x), lim z(t,x) = v(x) , 

t — >oo 

where the limit holds uniformly for x in bounded sets. To construct a(t,x), one can first define a 
smooth function a = d(t, x, z) such that 

ait -x z)-d(tx z)-{ 1 if l^-^)l> 2 A, 
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Then we solve the quasilinear Cauchy problem 



Zf — ct{t, x, z)z xx , 



z(0, x) = w(x) 



and set a(t, x) = a(t, x, z(t, x)) . We now claim that 

d , ,00 rx+p 

dt \J -00 J x-p 
Indeed, calling <j> = z x > and using (2.19) we compute 



z x {x) z x (y) dydx < 0. 



(2.20) 



d_ 
~dt 



CO fX-\-p 



/OO P 
-co J X 



(x) (j)(y) dydx 



co •/ x — p 

CO 



CO 

00 



t = (a(t,x)(f> x ) x , 

CO px-\-p 



CO J x — p 



(a 4> x (x)) <j>{y) + 4>(x) (a 4> x (y)) dydx 



/CO 
[a 4> x (x + p) - a<t> x (x-p)] (j)(x) dx 
-00 

/OO 
<f>(x) [a(j) x (x + p) - a<f> x (x - p)] dx 
-00 

acj) x {x) [<j)(x - p) - 4>{x + p)] dx 

-00 

<0, 

because 4>(t, •) is an even function, non-increasing for x > 0. From (2.20) it follows 

/oo fX+p poo fX+p 

I v x (x) v x (y) dydx < / / w x (x) w x (y) dydx . 
-ooJx — p J — 00 Jx — p 



□ 



Lemma 3. Let u be a solution of (1.1) defined for t 6 [0,r] and let w = w(t,x) as in (2.9)-(2.10). 
Set 

(2.21) 



a = -fil + (M) + «[Q(«(0)) - Q(u(t))] 



and let 



v(t, x) 



x/t if \x\/t < a , 

sgn(x) - a if \x\/t > a , 



(2.22) 



be a solution of Burgers' equation consisting of one single centered rarefaction wave of strength 2a. 
Then 

PT P CO PT P CO 

/ / [w(t, x+p)—w(t, x—p)]w x (t,x)dxdt<2 / / [v(t , x+p)—v(t, x—p)]v x (t,x)dxdt. 

JO J -00 JO J -00 

(2.23) 
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Proof. To compare the integrals in (2.23) a change of variables will be useful. We define (fig. 3) 



£€[0,<t], te[0,r\. 



For t G [0, r] and Q(t) — Q{t) < £ < a, we also consider the point y(t,£) > implicitly defined by 



w 



(t,oo) -w(t, y(i,£)) =(7-£. 



Notice that y(t,£) is defined only for t G , r] , or equivalently £ G [£(£) , <r] , where 

at) = K[Q(t)-Q(r)], t(0 = inf{t>0; [Q(i) - Q(r)] < ^} . 

For < £i < £2 < ^ and s > we have 

>(6-6)^ 

= ar(s,£ 2 ) -x(s,£i). 
Observe that, since w is odd and non-decreasing, 



(2.24) 



w 



+ 



(*, V ~ P) = max {w(t, y - p) , 0} = to(i , max{y - p , 0}) . 



Of course, the same is also true for v. Calling I w , I v the two integrals in (2.23) and using (2.24) 
at the key step, we obtain 



T PC 



T — 2 




/o Jm 

f'T f'(T 



d£dt 



Jtfc) 



■/ ^h/(t,6)-y(t,6) \<p 



w(t, y(t,0+p) -w(t, y(t,0~p) 
c-l/ / w(t, y(i,£)+p) y(t,0~p)\ d^dt 

d£id& dt 

= 4 J J meas {iG [0,r]; |y(i,£i) - y(t,6)| < p} d£id£ 2 
< 4 //meas{ f G [0,r] ; |ar(i,£i) - x(t,£ 2 )| < p} ^6^2 
>(t, x(t,0 + p) ar(t,0-p)l ded* 




jo 



< 2I V . 



□ 
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figure 3 



Corollary 1. Assume that all characteristic fields for the system (1.1) are genuinely nonlinear. 
Let u be a solution with initial data u(0, x) = u(x) having small total variation. Then, for every 
t, 5 > 0, the measures n\ + of positive waves in u(t,-) satisfy the estimate 

n ,. T 

\x-y\ <<5})dt = 0(l)-(ln(2 + r) + |lnJ|)<5-Tot.V a r.{n}. (2.25) 

i=i -* 

Proof. By Proposition 1 and the previous comparison lemmas, for every i = 1, . . . , n the integral 
on the left hand side of (2.25) has the same order of magnitude as in the case of a solution with a 
single centered rarefaction wave, of magnitude a = Tot.Var.{-u} < 1. Looking back at Example 1, 
from (2.3) we thus obtain 

/ Vj + ® tf) ({ (X, y); \x - y\ < 25}) = 0(1) ■ 5a (l + In ^) 

= 0(1) • (ln(2 + t) + \]h6\^6- Tot.Var.{n} . 

□ 

Remark 3. All of the above estimates refer to an exact solution u of (1.1). If u v — > u is a 
convergent sequence of front tracking approximations, the corresponding measures of i-waves in 
u v (t, •) converge weakly: fj, l v t — ^ \i\ for all % = 1, . . . , n and t > 0. Unfortunately, this does not 
guarantee the weak convergence of the signed measures 

For example (fig. 7), on a fixed interval [a, b] every u u might contain an alternating sequence of 
small positive and negative waves, that cancel only in the limit as v — > oo. However, by a small 



14 



modification of these front tracking solutions u v one can achieve the weak convergence (2.27) for 
each t in a discrete set of times {jr/N ; j = 0, 1, . . . , N}, with iV » e _1 . As a result, we obtain 
an arbitrarily accurate front tracking approximation (still called u) satisfying an estimate entirely 
analogous to (2.25), namely 




a a op\ )dt = 0(l)- ln(2 + r) + |ln£| v^-Tot.Var.{u}, (2.28) 



where we replace the S in (2.25) by for the application in Section 4. Here TZi denotes the set 
of rarefaction fronts of the i-th family and summation is over all possible pairs, including the case 
where the two indices a, (3 coincide. 



3 - Estimates on shock fronts 

We begin by estimating the sum in (1.9). The approximation v is discontinuous precisely at 
those times tj where an interaction occurs involving a large shock. Indeed, at such times the left 
and right states u~ , it+ across a large shock located at x = x a suddendly change. As a consequence, 
the viscous shock profile connecting these two states is modified. The two smooth functions v(ti—) 
and v(ti) will thus be different over the interval [x a — \fe, x a + y/e\ . To estimate the L 1 norm of this 
difference, the following elementary observation is useful. Given a smooth function (p = 4 > ( cr i <7 ')i 
its size satisfies the bounds: 

if 0(<t,O) = O for all a, then <j>(a, a') = 0(1) ■ \a'\ , 

if (j)( a , 0) = 0(0, a') = for all a, a', then <p(a, a') = 0(1) -\aa'\. 
We now distinguish various cases. 

1. At time ti a new large shock is created, say of strength \a a \ > p/2. In this case, since the new 
viscous shock profile is inserted on an interval of length 2-y/i, we have 

\\v(ti) -viU-)^ =0(1)-Vi\a a \. 

According to our construction, every large shock not present at time t = must grow from a 
strength < p/2 up to a strength > p at some later time r. Therefore, the sum of the strengths of 
all large shocks, at the time when ti when they are created, is 0(1) • So, where So = Tot.Var.ju}. 
The total contribution due to these terms is thus 0(1) • y^sSo- 

2. At time U a large shock is terminated. Since every large shock must have strength > p at some 
time and is terminated when its strength becomes < p/2, every such case involves an amount of 



15 



interaction and cancellation > p/2. Therefore, the total contribution of these terms to the sum in 
(1.9) is again 0(1) • ^e5 Q . 

3. A front 073 of a different family crosses one large shock a a . In this case we have 

\\v(U) - v(ti-)\\^ = 0(1) -y/i\a a \ \a p \ . 

These terms are thus controlled by the decrease in the interaction potential Q(u). Their total sum 
is 0(1)-^. 

4. A small front 073 of the same family impinges on the large shock a a . In this case we have 

H*i)-«(ti-)|| L i =0{l)- y/i\<rp\, 

Since any small front can join at most one large shock of the same family, the total contribution 
of these terms is 0(1) ■ \fe Sq. 

5. Two large /c-shocks of the same family, say of strengths a a , erg, merge together. In this case 

\\v(ti) - u(tj-)|| Ll = 0(1) • min{|cr a |, (cj^I} . 

As will be shown in (3.23), all these interactions are controlled by the decrease in a suitable 
functional by noticing that \a a \, \ap\ > 2y/e\ lne\. The sum of all these terms is thus found 

to be 0(1) • 5 Ve\lne\ . 

Putting together all these five cases, one obtains the bound (1.9). 

Next, we need to estimate the running error in (1.8) related to the big shocks, namely 

/•t [■x a + \/e 

Ebs = \v t + A(v)v x - ev xx \ dxdt . (3.1) 

Here the summation ranges over all big shocks in v(t, •). 
We first consider the simplest case, where the interval 

I a (t) = [x a (t) - 2^1, x a (t) + 2v^] (3.2) 

does not contain any other wave-front. In this case, observing that 

[A(ul(s)) - A«] - £ = ' 
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and recalling (1.19)-(1.20), the error relative to the shock at be written as 



E a (t) 



— v/e/2 ny/e N 
+ / 



d 



ds 2 



(3.3) 



Using the bounds 
d 



8s 



0(1) 



\sa a \/e 



ds 2 



(3.4) 



from (1.20) we deduce 

E a (t) = 0(l). |^exp|-^0(o}-(^0'(O + ^0"(o) dt 



'Ve/2 
esfl 

0(1) • / exp 

Jy/i/2 



\<Jc. 



0"n 



0(1) • / exp 



0(1) > a | exp{- 



|0Vv I 5 I . . o q ds 



+ 



2k„ 



2y^ J ^ ' / ' 

Since by assumption |cr a | > p/2 = c l\fe \ lne|, the above estimate implies 

£«(*) = 0(1) -e(l + |lne|) \a a \. 



(3.5) 



In the general case, our error estimate must also take into account the presence of other wave-fronts 
within the intervals I a (t). Indeed, for every point x a where large shock is located, we have 



E, 



m = I 

J X 



\v t + A(v)v x — ev xx \ dxdt 



x a —\fe 

0(l)- £ (l + |ln £ |) \a a \ +0(1) 



eI a (t),\xf,-Xy\<2y/Z x e &I a (t),e€BS 



In the following, we introduce three different functionals, which account for: 

• products IcTq.0"^! of fronts of different families, 

• products (ctqit^I where a a is a large shock and 073 is a rarefaction of the same family, 

• products | (TaiT/? | of shocks the same family. 

By combining these three, we form a functional Q(u) such that the map t t-> Q(u(t)) is 
non-increasing except at times where a new large shock is introduced. Moreover, the total in- 
crease in this functional at times where large shocks are created will be shown to be 0(1) • 
y/e\ lne|Tot.Var.{n}. 
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We begin by defining 

where the sum extends over all couples of fronts of different families (small shocks, big shocks, 
rarefactions). The weights W^p G [0, 1] are defined as follows. If hp < k a , then 



if xp < x a — 2y/e, 

if xp G [x a — 2yfi , 
if xp > x a + 2y/e. 



W U = \ \ + X \jT if xpe[x a -2^I, x a + 2^/I\, 



If instead kp > k a , we set 



1 if xp < x a — 2y/e, 



Wlp=l \ - if xp€ [x a -2Vi, x a + 2Vi\, 

if xp > x a + 2\fe . 

By strict hyperbolicity, we expect that the functional Q b will be decreasing in time. Indeed, its 
rate of decrease dominates the sum 

Y WaVpl, 

k a ^k(}, \x a -xp\<2y/e 

containing products of nearby waves of different families. 




figure 4 

Next, given a big shock a a of the fc Q -th family located at x a , we write: 
1Z a to denote the set of all rarefaction fronts of the same family k a , 
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S a to denote the set of all shock fronts of the same family k, 



To control the interaction between large shocks and rarefactions of the same family, we define 
the weight 



and the function (fig. 4) 



w a (x) = < 



^(x)=min^ + ^-J 



E 

P&7l a , x p e[x, x a ] 

E °p 



(3.8) 



if x > x c 



Calling 



u) Q (x) 



f-|<7 a |/4 if w a (x) <-\a a \/4, 

w a (x) if |u> a (x)| < |<r a |/4, 

cr a |/4 if w Q (x) > |cr a |/4, 



we then define 



o"(«)= E / w£(*)A 



(3.9) 



a€BS 



By using the function with cut-off w a , instead of w a , in (3.9) we are taking into account only the 
rarefaction fronts ap of the same family k a , such that the total amount of rarefactions inside the 
interval [x a , xp] is < |<r Q |/4. If no other fronts of different families are present, this guarantees that 
all these rarefactions op are strictly approaching the big shock a a . Indeed, the difference in speed 
is \xp — x a \ > \a a \/4. As a result, the functional Q\u) will be strictly decreasing. On the other 
hand, if the interval [x a ,x^] also contains waves of different families, the above estimate may fail. 
In this case, however, the decrease in the functional Q^(u) compensates the possible increase in 
Q\u). 

Finally, to control the interactions among shocks of the same family, for each shock front a a 
(of any size, big or small) located at x a , we begin by defining (fig. 5) 



z a (x) 



E 



E 



3ov 



/3e<s a , x<xp<x a 



P£TZ a , x<xp<x a 



+ 



E 



E 



3a r - 



if OC *\ Q; } 

if x > x a . 



0eSn, X a <Xp<X 



fielZc,, X a <Xfj<X 



Then we set 



z a (x) 



min {z a (x') ; x < x' < x Q } if x < x a , 

^ max [z a (x') ; x a < x' < x} if x > x a . 
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figure 5 



Notice that z a is a non-decreasing, piecewise constant function, with (x — x a ) z a (x) > for x ^ x a . 
Using the weights 

f r ~ i \i _1 




if x < x c 
-1 . r 

if X > x c 



we now define 



Q\u) = Kl / W^(ar) W^A^ . (3.10) 

Notice that in this case the summation runs over all shock fronts. If ap is a shock located at 
X/3, then [w!(x/3)] 1 roughly describes the amount of shock waves inside the interval in 
excess of three times the amount of rarefactions. If the interval [x a , xp] does not contain waves of 
other families and the function x i-> z a (x) has a jump at at x = xp, then the two shocks a a , ap 
are strictly approaching, hence the functional Q"(tt) will decrease. On the other hand, if waves of 
different families are present, the above estimate may fail. In this case, however, the decrease in 
the functional Q b (u) compensates the possible increase in Q"(«). 

In the definition of z a , notice that the strength of rarefactions is multiplied by 3, to make 
sure that couples of shocks a a , ap entering the definition of Q^(u) are always approaching each 
other (except for the presence of fronts of different families in between). An example is shown in 
fig. 6, where two nearby shocks move apart from each other because there are sufficiently many 
rarefaction waves in the middle. Because of the factor 3, the function x z a (x) will be constant 
at the point xp. Hence the product IfJaCr^l will not appear within the definition of Q^(u). 

We now consider the composite functional 

Q(u) = Vi| hie| ■ (Ci T(u) + C 2 Q\u) + C 3 Q*(u)) + v^Q S («) • (3.11) 
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x a b x 

figure 6 figure 7 



Here 

T(u) = V{u) + C Q{u) (3.12) 

is a quantity which is decreasing at every interaction time. Its decrease dominates both the amount 
of interaction and of cancellation in the front tracking solution u. Observe that 

Q{u) = 0(1) • y/e\ hie | Tot.Var.{u} . (3.13) 

Indeed, by the definition of W|(x), we have 

r ,.Tot.Var.{«} -, 

W*(x)D x z a = 0(1) ■ ——ds = 0(l)- \lne\. (3.14) 
J Jo s + e 

Using (3.14), it is now clear that 

Q\u) = Oil) ■ | lne| Tot.Var.{n} , T(n), Q b {u), Q\u) = 0(1) • Tot.Var.{u} . (3.15) 
The bound on (3.11) now follows from (3.15). 

Lemma 5. For a suitable choice of the constants C\ » C2 » C3 » 1, if Tot.Var.{u} remains 
small, then at each time t* where an interaction occurs the following holds. If a new large shock 
of strength \a a \ > 2^/e\lne\ is created, then 

AQ = Q(t+) - Q(r-) = 0(1) • ln £ |K| . (3.16) 

If no large shock is created, then 

AQ<0. (3.17) 

Proof. Notice that the weight a is always < 1. For a newly created large shock a a , the 
increase in the functional Q b (u) can be estimated as 

AQ b (u) = 0(1) • |<7 a | Tot.Var.{u}. (3.18) 
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Similarly, since W\ < 1, it is clear that the increase of Q^(u) due to a new large shock a a is 

AQ\u) = 0(l)-\a a \. (3.19) 
The estimate on the increase of the functional Q}{u) is different. In this case, the integral 

W*(x)W*(x)D x z a 



I 



is bounded by 



Hence, 



,.Tot.Var.{u} 
Oil ) - I __ dx = 0(1)| lne|. 



o 



AQ s (u) = Oil) ■ \a a \ |lne|. (3.20) 



Together, (3.18)-(3.20) imply (3.16). 

Next, we prove (3.17). Assume that at time t* an interaction occurs without the introduction 
of any new large shock. We will show that the functional Q(u(t)) decreases. 

First we look at the change in Q^(u) and Q\u). Since the weights and are uniformly 
bounded, it is straightforward to check that the change in these two functionals at time t* is 
bounded by a constant times the decrease in the Glimm functional T(u(t)) in (3.12). Hence, by 
choosing C\ » C2 » C3, the quantity 

CiT (u) + C 2 Q ] '(u) + C 3 Q\u) 

is not increasing in time. 

The analysis of Q^(u) is a bit harder. We will show that the change of Q^(u) at the interaction 
time t* is of the same order of magnitude as | lne|| AT(«)|. Here and in the following, AT denotes 
the change in T(u(t)) across the interaction time. As a preliminary, we notice a basic property of 
the weight function W|(x). For any fixed location x = xq, we have 

J2K\Wl(x ) = O(l)-\\n£\, (3.21) 

aes 

E Kl / (Wt(x)) 2 D x ~z a + Yl Kl / {Wt{x)) 2 D x ~z a 

aGS,x(a)<x J x ° a€S,x(a)>x J 00 (3.22) 

= 0(1) ■ \\ne\. 

The proof of the estimates in (3.21) and (3.22) is straightforward by noticing that the functions 
fi x ) = i+7> d( x ) = (x+e) 2 are convex anc l bounded away from zero for x > 0. And the left hand 
sides of (3.21) and (3.22) are bounded by the following single and double integrals respectively: 

/•Tot.Var.{«} /.Tot.Var.{«} ^Tot.Var.M 

/ f(x)=0(l)\\ne\, / / g(y)dydx = 0(1) ■ \\ne\. 

JO JO Jx 
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Now we are ready to estimate the change in Q$(u(t)) at time t*. Note that, in some cases, it 
is possible that the interaction does not change the functional. In the following, we will consider 
the case where Q"(«) does change across the interaction. Depending on the types and families of 
the waves involved in the interaction, we have the following four cases. 

1. Two shocks of the same family interact. Let f3\ and /?2 be the two interacting shocks, say of 
the i-th family, and call (3 the outgoing i-shock. We also let ol\ and «2 be any two shock fronts on 
the left and right of the interaction point respectively, so that x ai < xp < x a2 at time t = t* . For 
any shock front a a at time t* , set 



Qi = Wa\ J W*(x)W*(x) D x z 



Observe that 



Q\ - « + <4) < - i xriT + • 1 lne| |AT| • (3 - 23) 



Indeed, 

<7 / 3=<7 /3l +<7 i 9 9 +0(l)-|AT|. 

Moreover, recalling (3.8), we see that after the interaction we lose the term 



™k (*& ) {Wl ( Xf s 2 ) + Wl (x Pl ))\ a Pl ap, 



> 



l°"/3il + l°73 2 l + e 

Notice that W^(x) (respectively W| (cc)) does not change across the interaction for x < xp 
(x > xp). The change in the Q^ ai , i = 1, 2, can be estimated as follows. When «j is of the same 
family of /3j, = 1,2), by (3.22) we have 

AQ» =( Igglf^J )wHx(8)) 

1 \e + \a ai \ + \a p \+I e + \a ai | + \a h | + / £ + 1^1 + 1^1 + 1^1+// " 

/» CO 

+ 0(1)-|AT||<t Qi | / (^J 2 (x)^z Ql 

POO 

<0(l)-|a Ql ||AT|^l i (x /3 ) + 0(l).|AT||a Ql | / {Wlf(x)D x ~z ai . 

Here and in the following, we assume that the whole strength of i = 1, 2 and of /3 appear in the 
functional . Moreover, / represents the sum of the strengths of the i-shocks between f3\ and 
cci that appear in Q^ ai - The other cases when part or none of the above wave stengths appears in 
Q^ ai can be treated similarly. 

By summing over ot\ and using (3.21) and (3.22), we find that the total change of Q^ ai is 
0(1) • | lne||AT|. A similar estimate holds for a^- 
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Now consider two shock fronts a±, «2 of the j-th family, with j ^ i. Notice that the change 
of the weight function W|. (x), i = 1, 2, is at most of the order of (W|. (x)) 2 | AT| when x lies on 
the opposite side of x ai w.r.t. xp. Together with (3.22) this yields 

X>T||<t Qi | / {wl i )\x)D x ~z ai +Y J \^\W a ,\ / [Wl)\x)D x ~z a2 
= 0(1) ■ |lne||AT| . 

If 7 is a newly created shock of the j-th family, then the new term has size 0(1) ■ |cr 7 | | In e\ . 
Hence, the total sum of these new terms over 7 is of 0(1)| lne||AY|. And this completes the 
discussion on this case. 

2. Interaction of a shock with a rarefaction front of the same family. Let [3i and j3 2 be a shock 
and a rarefaction front of the i-th family, interacting at time t* . 

First, consider the case where the shock (3\ is completely cancelled and hence the decrease in 
T(u) is of the same order as j3\. In this case the term disappears after the interaction. Let a.\ 
and «2 be shock waves of the j-th family on the left and right of the location of interaction. For 
both cases when i = j and i 7^ j, by (3.22) we have 

X>T||<t Qi | / (I^J^A^+^TlATII^I / (Wl 2 ) 2 (x)D x z a2 

ai Jx P a 2 J -°° 

= 0(1) • |lne||AT|. 

The same argument applies to the change in Q^, related to the newly created shock 7 of the j-th 
family, when j 7^ i. In this case, the total change in is again 0(1) • | lne||AT|. 

In the case where the interaction produces an outgoing i-shock so that the rarefaction /3 2 
is completely cancelled, the analysis is as follows. First, notice that the increase in Q"(«) due to 
the newly created waves is 0(1) • | lne||AT|, with |AT| = 0(1) • \&p 2 \. 

Next, the difference between and comes from the changes in (W^) -1 and zp 1 which 
are at most of the order of ap 2 at each x. Hence 

Qh-Qh < j wl(x)wl(x)D x ~z^ - j wl{x)Wl{x)D x ~zp,) 

/•Tot.Var.{«} , ,| CT/32 | , 

= 0{l)-\ap 2 \\ap 1 \ - j-^j +0(l)-\a h \ ^ 

Pl Jo (e + k/3j+y) 2 01 Jo e + la^l+y 

= 0(1)> & | = 0(1).|AT|. 

The change in Q^ a . also comes from the change in W|. (x) and z ai (x). Since the weight W|. (x) 
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decreases as x moves away from x ai , we have 

Ewkj/ (wm X )D x z ai +j2mk 2 \ / (wif(x)D x 

+ 0(l)-\a^(^\a ai \Wi 1 (x fh ) + Y / W a2 \Wi 2 (x fjl ) 



«i 



= 0(1) • |lne||AT|. 

3. Interaction of a shock and a rarefaction front of different families. To fix the ideas, let (3\ be a 
shock of the z-th family and /?2 be a rarefaction wave of the j-th family with i > j. Assume /?i and 
P 2 interact at time t* and denote the outgoing wave of the i-th family by fii, and the j-th. family 
wave /?2- Moreover, let 7 be a newly created shock front of the k-th family, k 7^ i, j. By a standard 
interaction estimate, we have 

K-^J=0(1)-|AT|, |<r 7 | = 0(l)-|AT|, i = l,2. 

Thus, if we consider two shock waves c^, i = 1, 2 of the A;-th family located on the left and right of 
the interaction point respectively, as in the analysis of Case 2 we have 

(/■OO 
X>T||<t Qi | / {Wl) 2 {x)D x ~z ai +Y,\^\W a .\ / (Wl 2 ) 2 (x)D x ~z a2 

= O(l) • |lne||AT|. 

In addition, 

\Ql - Qjj = 0(1) • I lne||AT|, |Q^| = 0(1) • | lne||AT|. 
Here we assume that j3\ is a shock wave. In the other case, we have Qt = 0(1) • | lne||AT|. 

4. Interaction of rarefaction fronts of different families. The change of Q^(u) in this case only 
comes from the new fronts created by the interaction. Therefore, as in the analysis of Case 3, the 
total change in Q$ is bounded by 0(1) • | lne||AT|. 

Based on the analysis of the above four cases, we see that by choosing C\ to be sufficiently 
large, then the nonlinear functional Q(u) is non-increasing at the interaction time when no new 
large shocks are introduced. This completes the proof of the lemma. D 



4 - Proof of the main theorem 
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Relying on the analysis of the two previous sections, we can now conclude the proof of The- 
orem 1. We briefly recall the main argument. If one defines the mollification v 5 = u * cps with 
5 = y/e, the estimates (1.7) hold, while (1.13)-(1.14) imply 



J J\v s t+ A(v s )v s x -ev s xx \dxdt = 0(l)- J T Josc.{u; [y - 5, y + 8]} \du(y)\ 

r ° ( 4 -!) 

= o(i)- / y, \°<*°p\ dt 

\x a {t)- Xfi {t)\<5 



dt 



In this case, the presence of big shocks gives a large contribution to the right hand side (4.1), 
namely 

f £ \a a (t)\ 2 dt (4.2) 
Jo aeBS 

To get a more accurate estimate, in a neighborhood of each big shock we replaced the mollification 
with a (modified) viscous travelling wave, according to (1.21). By doing this, we picked up more 
error terms, namely: 

• The terms related to the interactions of big shocks with other fronts. The analysis at the 
beginning of Section 3 has shown that the total contribution of all these terms satisfies the 
bound (1.9). 

• The errors due to the difference between the rescaled profiles u> a in (1.19) and the exact 
travelling wave profiles u a . According to (3.4), the total strength of these terms is 

f Yl E a (t) dt = O (1) ■ re + \ In e\) Tot. V&r.{u}. (4.3) 

On the other hand, we removed the contributions of all terms in (4.2). For the function v 
defined at (1.21) we thus have 



J J \v t + A(v)v x — ev xx \ dxdt = 0(1) • re(l + | lne|) Tot.Var.ju} 

+o(i). /Y j2 i<wi- E K\ 2 )dt. 



(4.4) 

\cr a \ z I dt . 

■\x l3 -x 1 \<2^/E a€BS 

The main goal of this section is to show that the last integral in (4.4) can be estimated as 



f T ( Yl i<wi - E i^i 2 ) dt = ■ E f ( E i^i) dt 

V \x p -x-,\<2^ aeBS ' i=l Jo ^ /3,-yeni,\xi3-x^\<8^i ' (45) 



+ o(i)- f 

Jo 



j t Q«t)) 
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dt + 0(l) ■ v£|]ne|rTDt.Var.{u}. 



Using the estimate (2.28) on the spreading of positive wave-fronts and the bounds (3.15)-(3.17) 
concerning Q(u), from (4.5) we obtain 

j ( Y \ a P a i\ - Y l "' 2 ) dt = ' ^ + r )v / e|lne| • Tot.Var.{u} . 

^ \xp-x. t \<2^ a&BS ' 

This will complete the proof of the estimate (1.3). 

The remaining part of this section is devoted to a proof of (4.5) which is a consequence of the 
following lemma. 

Lemma 6. Outside interaction times, one has ^Q(u(t)) < and 

Y w^i - Y I ct «i 2 = °( 1 )-Y ( Y h^i) 

\x p -x^\<2^ a€BS i = l ^ ,\x l3 -x- l \<8^ ' (45) 



+ 0(1)- 



l«(«(t)) 



+ C?(l)-Ve| lne|Tot.Var.{n} 



To help the reader work his way through the technicalities of the proof, we first describe the 
heart of the matter in plain words. 

After removing the terms in (4.2) related to large shocks, the left hand side of (4.6) still 
contains the sum 

Y i a °i 2 ' 

a£<S<S 

where SS denotes the set of all small shocks. According to (1.22), the maximum strength small 
shock is < 4y/e | In e| . Hence the above sum is estimated by 0(1) • \fe \ In e| Tot.Var.ju} . 

Next, consider any interval J of length We first estimate the restriction of (4.6) to fronts 
inside J, i.e. 

@= Yj ka^l — Yj \ a a\ 2 ■ 

x a eJ, \x a -xp\<2y/e x a £J,aeBS 

It is convenient to split 6 into various sums: 

© b = Y Wa<7p\, @i = Y WaO-p\, 

\<2y/e, k a ^kf3 \<iy/l, aeSi, pent 

©r ef = Y i^i, ef= y K^i- 

x a €J, \x a —xp\<2y/e, a,/3e7li x a eJ, \x a —xp\<2y/e, a,/3eSi 

If b dominates all other terms, then the whole sum can be controlled by the rate of decrease in 
the functional Q b , related to products of fronts of different families. The alternative case is when J 
contains almost only waves of one single family, say of the i-th family. If ®\ is the dominant term, 
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then is controlled by the decrease of the functional and QK If 6| dominates, then J contains 
mainly z-shocks, and is controlled by the decrease in QK Finally, if 6J aref dominates, then there 
is nothing to prove, because the sum over all couples of nearby rarefactions appears explicitly also 
on the right hand side of (4.6). 

Covering the real line with countably many intervals Jg of fixed length, we eventually obtain 
the desired result. 

Proof of Lemma 6. Since the system is strictly hyperbolic, the definition of the functional 
Q b (u) implies 

(4.7) 

\ l/Trt/T I 



\xf3—x~ f \<2^/e, kp^k 1 

for some constant C2 > related to the minimum gap between different characteristic speeds. 
Hence the terms containing a product of two waves of different families on the left hand side of 
(4.6) are controlled by the decreasing rate of Q^(u). 

In the following, we only need to show that the products involving one or two shock waves of 
the same family can be controlled by the decreasing rate of the nonlinear functional Q(u), plus the 
quantity in (2.28) and y/e\ lne|Tot.Var.{-u}. 

By the definition of Q^(u), we know that the rarefaction waves located in I a (t) involved in 
Q^ a (u) approach to the large shock wave a unless there are waves of the other families in between. 
Hence, if we use BS' to denote the set of big shocks a such that the total strength of small wave 
fronts within the interval I a (t) is < ^j^, then for a € BS' , we have 

<r C3 ^ I lx °M I, 

< \<J a <J(3\ + —rT WaOph 

v x ei a (t),i3en a v p,x p ei a {t),k p ^k a 

where 

^ r Wr 



Qi(u) = J W*(x)D a 



On the other hand, the functional Q^{u) is defined for all shock waves no matter they are 
small or large. In this way, its time derivative yields mainly the product of two shock waves of the 
same family with distance < 2-y/e. Let 



Qi(u) = \a a \ J W*(x)W*(x)D a 



Since there is a factor 3 in front of the summation of rarefaction waves in the definition of Q* (u) , 
this guarantees that all the shock waves appearing in Q$ a approach to the shock wave a if there is 
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no waves of other families in between. Notice that there is a constant e in the denominator of the 
weight function W^(x). Thus, for any shock a G S, 



d_ 

dt 



c 5 



Ql(u) < % V . \a a a \ + .... 



E 



f)en a ,x p ei a {t) 



0(1) 



E 



^-■t= E 



+ 



0(1) 



E 



pen (<) 
\a a ap\ + C(l)A/i|<r a |, 



(4.9) 

where c 4 ,c 5 > are constants independent of e, I a (t) = [x a — 2y/e, x a [\J ]x a , x a + 2y/e\, and J2 
means that the summation is over all shocks (3 with the property that the total strength of all 
shock fronts between a and (5 with xp G I a {t) is > e. 

By noticing that the time derivative of T(u) is zero outside interaction times and by choosing 
Ci » Cs » 1, based on the estimates (4.7)-(4.9), the increase of (4.6) can be given as follows, 
by considering separately the products involving large shocks, and those involving only small wave 
fronts. 

For a large shock front a, consider the summation 



E l^a^l- 

Since the sum of all products lowers) when k a ^ kp is controlled by (4.7), we have 



(4.10) 



(4.11) 



(4.12) 



E < 0(1) ■ 

xpei a (t) 

provided that waves of different families dominate, say 

E \°e\>\ E M- 

It thus remains to consider the case when (4.12) does not hold. We then have 

E m - 4 E M- 

xp£l a (t),kp=k ct (3,x g el a {t) 

In this case, if a G BS 1 , then the summation of lo-aO^I for j3 G 1Z a U S a is controlled by (4.8) and 
(4.9) together with (4.7). Therefore 



(4.13) 



, n 

E E 1^1 < ^(i)- (E E 

*eBS' xg ei a (t) i=1 , \x p -x^\<i^l 



(4.14) 



+ 0(1) 



+ C(1) -eTot.Var.{u}. 
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Moreover, by (4.9), for a G S, if 

£ E W ( 4 - 15 ) 

pes a ,xpei a (t) (3eK a ,Xf)ei a {i) 

then 

£ < C?(l) ■ ^- ^(CaVillnelQS, + ViQg.) - (4-16) 

For a large shock wave, it now remains to consider the case when a 6 BS" = BS — BS' 
satisfying (4.13) and 

E E M- 

pes a ,x p ei a (t) (3eK a ,Xf)ei a {t) 
We denote the set consisting all these large shock waves by BS'" . Notice that this is a subset 
of BS" . Roughly speaking, for a 6 £>£'", the small wave fronts is not small compared to a and 
rarefaction waves of fc a -th family dominate in I a (t). Hence, for a € BS'" , one has 

E E i<wi<o(i)- E i^i, (4-i7) 

which is controlled by the corresponding part of (2.28) in the interval [x a — 4-y/i, x a +4^/e\. Hence 

E E We.*?] < o(i) ■ e i^i' ( 4 - 18 ) 

which is estimated by (2.28). 

Combining (4.11), (4.14), (4.16) and (4.18) we obtain 

E E km<o(i)-(e E 

a€BS Xf3 ei a {t) V * =1 PneKi ,\ Xp -x^\<&y/I 7 (4.19) 



+ 0(1 > * 

Now it remains to show the sum of products of small wave fronts of the same family satisfies 
the same bound: 

E k^i<o(i)-(e E i^i 

x,PeSSun,\x a - X0 \<2^,k a =k ^ i=l P,-y^i , \xf,-x y \<8y/e (4.20) 



j t Q{u{t)) +0(1) -eTot.Var.{n} 



+ —Q(u(t)) + v^|lne|Tot.Var.{u} 

To obtain the estimate (4.20), we divide the real line into a union of closed intervals of length 
2y/e, i.e. 1R = [j i Ji with Jj = [2i^/e, 2(i + l)y/e] • We denote by and rf respectively the total 
strengths of small /c-shock and /e-rarefaction fronts contained in the interval in Jj. We have 

E Wa<Tp\<^2 E E WaVpl- 

a,PeSSUK,\x a -xp\<2,/e,k a =kp i aeSSUTZ,x a eJi peSSUTZ,\x a -x l3 \<2^e,k a =k f3 



30 



To estimate the quantity 



E E \wp\> ( 4 - 21 ^ 

a€SSl)K,x a €Ji /3eSSun,\Xc,-x \<2^/e,k a =k 



we consider the following two cases. 

• For a given k, > ^-(rf_ 1 + r\ + rf +1 ). In this case, from (4.9) we deduce 

I E q'm<-^ £ E E k^i 

aeSSk, x a € Ji a€SS k ,x a eJi p € s a , xpel a (t) 

+ ^ E E K^i+o(i)Vi E 

V aeSS k ,x a eJi Xpel a (t),kf3^k a a€SS k ,x a e.Ji 



aess k ,x a eJi pes 
0(1) 



(4.22) 



E K! 2 + 0(l)Vi E i^i 

V aeSS k ,x a €.Ji a€SS k ,x a £Ji 

E E |(7Q ^ | + ^7r S E 

+ 0(1)1 lne| E k„| + 0(l)V£ E 1^1- 

a£LSS k ,x a GJi a&SS k ,x a £Ji 

Here we have used the fact that cr a < 4-y/ij lne| for a € <S<Sfc. By (4.22) we see that those terms 
containing a product with a 6 S<Sfc and /3 G iS<Sfc in (4.21) can be controlled by ^Q{u) up to 
an error of the order of y/e\ lne|Tot.Var.{lZ}. Since the total stength of all small fc-shocks in Jj 
dominates the total strength of all /c-raref act ions in Ujii-i Jji the products of a £ SSk and 
(3 G IZk, and the products of a 6 IZk with (3 € 55^ for xp G Jj in (4.21), are also controlled by 
■^Q(u) up to an error of the order of yf£\ lne|Tot.Var.{n}. Moreover, those products of a € IZk 
and (3 G 7£fc in (4.21) are controlled by the corresponding parts of (2.28) the interval Uj=;-i Ji- 

Hence, it remains to consider the product of a G IZk and (3 6 SSk with xp € Ji_i U Jj+i. To 
fix the ideas, we consider the case when a 6 TZk, (3 6 55^ with a;^ G Jj-i, i.e., 

E (4.23) 

a£7l k ,f3£SS k , x a GJi, xpG.Ji-i, \<2^f£ 

When 4_i < T t (r*-2+ f f-i+ T f)> (4.23) is controlled by (2.28) in the interval U$tl-2 J J- Otherwise 



E ka^l < ( E MJ 



2 



aeTZ k ,/3ess k ,x a eJi,xpeJi-i, \x a -xp\<2y/e x p,xpeJi-i,pess k 

which can be controlled as in (4.22), using (4.7) to control the A;-shock fronts in Jj_i. 
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• Now assume that sf < ^f(rf_! +rf + rf +1 ). In this case the total strength of all /c-rarefaction 
fronts in Ujii-i Jj dominates the total strength of /c-small shocks in Jj. As done previously, we 
only need to consider the case when a 6 SSk U IZk and (3 G SSk with xp G Ji_i U Jj + i in (4.21) 
because all the other terms can be controlled by (2.28) in the corresponding interval U^i-i Jj- 
For illustration, we discuss the following two terms, 

E i^i' ( 4 - 24 ) 

|<2v^ 

and 

E kaff/jl, (4.25) 

a,/3eSS k , x a eJi, xpeJi-i, \x a —xp\<2^fi 

respectively as follows. The other terms can be handled similarly. 

Concerning (4.24), when s\_ x < ^■(rf_ 2 +rf_ 1 +rf), it can be controlled by the corresponding 
terms in (2.28) in the interval (Jj=i-2 Jj- Otherwise, when sj ; _ 1 > ^-(r^_ 2 + r\_ x + rf), we have 

E ^ ( E m) ' 

a€K k ,l3€SS k ,x cl €Ji,Xf 3 EJi- 1 , \x a -xp\<2^ ^ /3eSS k ,x g eJi-i 

which can be estimated as in (4.22), using (4.7) to control the fc-shock fronts in the interval Jj_i. 

Concerning (4.25), when s^_ 1 > ^-(rf_ 2 + r\_ x + rf), then a similar argument as in (4.22) 
can be applied, using (4.7). 

Otherwise, 

E K^I<(E^) 2 , (4-26) 

a,f3eSS k ,x a e.J i ,xi 3 eJ i _ 1 , |x„-x /3 |<2 v /e j—i-2 

which can be controlled by the corresponding term in (2.28) in the interval |J*i*_ 2 Jj- 

Notice that each interval Ji can be counted no more than three times. By combining (4.22)- 
(4.26), we have desired estimate on (4.21) for small wave fronts so that (4.20) holds. In summary, 
(4.19) and (4.20) imply (4.6), completing the proof of the lemma. 

□ 

Remark 4. In the proof of the error estimate (1.3), the three basic ingredients are: 

• The existence of uniformly Lipschitz semigroups of approximate (viscous) solutions. 

• The decay of positive waves, due to genuine nonlinearity, 

• The exponential rate of convergence to steady states, in the tails of travelling viscous shocks. 
Assuming that all characteristic fields are genuinely nonlinear, we thus conjecture that similar error 
estimates are valid also for the semidiscrete scheme considered in [Bi]. In the case of straight line 
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systems, based on the analysis in [BJ], it is reasonable to expect that analogous results should also 
hold for the Godunov scheme. 

Remark 5. In the case where all characteristic fields are linearly degenerate, solutions with 
Lipschitz continuous initial data having small total variation remain uniformly Lipschitz continuous 
for all times, as shown in [Bl]. Therefore, the easy error estimate (1.16) can be used. For systems 
having some linenearly degenerate and some genuinely nonlinear fields, we still conjecture that 
the error bound (1.3) is valid. A proof, however, will require some new techniques. Indeed, the 
contact discontinuities that may be generated by shock interactions at times t > can no longer 
be approximated by viscous travelling profiles. 
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